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Abstract 

We present analytic results for ground-state properties of Hubbard-type models in terms of the 
Gutzwiller variational wave function with non-zero values of the magnetization m. In dimension 
D = \ approximation-free evaluations are made possible by appropriate canonical transformations 
and an analysis of Umklapp processes. We calculate the double occupation and the momentum 
distribution, as well as its discontinuity at the Fermi surface, for arbitrary values of the interaction 
parameter g, density n, and magnetization m. These quantities determine the expectation value of 
the one-dimensional Hubbard Hamiltonian for any symmetric, monotonically increasing dispersion 
e^. In particular for nearest-neighbor hopping and densities away from half filling the Gutzwiller 
wave function is found to predict ferromagnetic behavior for sufficiently large interaction U . 
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I. INTRODUCTION 



Quantum-mechanical many-body problems can almost never be solved exactly. In this 
situation variational wave functions have proved to be particularly useful. Although they 
describe correlations among the particles only in an approximate way, they have the advan- 
tage of being explicit and physically intuitive. In particular, they allow for investigations 
even when standard perturbation theory is not applicable, or is untractable. 

Variational wave functions can, for example, be obtained by applying a suitably chosen 
correlation operator (e. g., the interaction part of the Hamiltonian under investigation) to 
a simple one-particle wave function. For the one-band Hubbard modelEJoH 

H = Y^ efcflL^fca + ^ 5Z ^^T^ii' (1) 

fccr i 

which is often used as an effective model to understand electronic correlation phenomena 
like itinerant ferromagnetism in transition metals, high-temperature superconductivity and 
the Mott-Hubbard metal-insulator transition, the simplest projected wave function is the 
Gutzwiller wave function (GWF)@ 

I^g) = g^^^'\%) = n - (1 - 3)d] |$o), (2) 

i 

where is a variational parameter (usually < (7 < 1), Z)j = fii^hii is the operator for 
double occupation at lattice site z, and the starting wave function |$o) is a product state of 
spin-up and spin- down Fermi seas 

l'^o)= n ^fcJvac). (3) 

fccr 

Using the GWF one may, in principle, calculate expectation values of any operator A as 
{A)q = (\1/g|v4|\1/g)/(\I/g|^g)- By the variational principle the energy expectation value 
{H)g is an upper bound for the true ground-state energy of H. 

The properties and quality of the GWF have been subject of detailed investigations (for 
an early review see Ref. A diagrammatic theory for the calculation of expectation values 
in terms of the GWF, valid in arbitrary dimensions D, was formulated by Metzner and 
Vollhardt.@ (Ref. ^ is hereafter referred to as MV.) In particular, for systems without net 
magnetic polarization m = n-|- — = (i. e., particle densities n^. = n/2, with fcpt = ^fj. 
in Eq. (^) they calculated the momentum distribution = {(''kcT^ka)G and the double 
occupation d = {J2i^i)G/L analytically in D = 1 for all values of g and n, where L is 
the number of lattice sites. The analytic calculation of correlation functions, in particular 
of the spin-spin correlation function, by Gebhard and VollhardtH showed that in the non- 
magnetic case, for U ^ 00, the results obtained with the GWF are in very good agreement 
with exact analytic and numerical results for the antiferromagnetic Heisenberg chain. In 
fact, HaldanJ and Shastryi discovered that the GWF at (? = is the exact ground state 
wave function of the antiferromagnetic Heisenberg chain with exchange. Results were 
also obtained in dimensions D = 1,2,3 using numerical techniquesi'0 and finite orders of 
perturbation theory.0'0 Within the diagrammatic approach of MV it also became possible 
to derive the well-known Gutzwiller approximation in the limit of infinite spatial dimensions 
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{D = oo).ll^ Comprehensive investigations in this hmit were made possible by the approach 
of Gebhard,0 which allows for explicit evaluations of expectation values for arbitrary start- 
ing wave functions |$o) (including ones with broken symmetry) and facilitates the expansion 
in 1/D around D = oo. This approach was also extended to multi-band Hubbard models;lll 
recently that method was combined with density functional theory, and applied to ferromag- 
netic transition metals.li^ The Gutzwiller approximation also describes a correlation-induced 
transition from metal to insulator, the Brinkman-Rice transition.EJ We recently investigated 
the effect of correlated hopping, whichfor the GWF can be calculated in terms of Uka and 
d in any dimension, on this transition.ta 

Up to now the analytic calculation of expectation values in D = 1 was limited to the 
unmagnetized paramagnetic phase. In view of the renewed interest in the microscopic foun- 



dations of metallic ferromagnetism (see Ref. 19 for a review), it is desirable to perform such 



evaluations also for the GWF with non-zero magnetization (m 7^ 0). In this paper we show 
that, in spite of formidable technical complications, it is indeed possible to evaluate such 
expectation values even for finite magnetization. From suitable canonical transformations 
we obtain diagrammatic relations and reduce the expectation values for m 7^ to those 
for m = 0. Thereby it becomes possible to calculate the double occupation d and the mo- 
mentum distribution Uka for arbitrary values of the correlation parameter g, density n, and 
magnetization m. Furthermore we derive an expression for Uka in closed form, which was 
not available up to now even for zero magnetization. These quantities determine the energy 
expectation value and thus the optimal variational parameter and spontaneous magnetiza- 
tion. 

The paper is structured as follows: In Sec. fl| the diagrammatic formulation is used to 
derive diagrammatic relations from canonical transformations, valid in arbitrary dimensions. 
The evaluation of expectation values is derived in Sec. The resulting magnetic phase 
diagram for the Hubbard model in D = 1 is presented in Sec. 0. The conclusion in Sec. [V] 
closes the presentation. 



II. DIAGRAMMATIC FORMULATION IN ARBITRARY DIMENSION D 
A. General formalism 

The double occupation d{g,n,m) and the momentum distribution nkcr{g,n,m) of the 
GWF are required for the calculation of the variational energy, Eq = {H)q/L. Another 
quantity of interest is the discontinuity of Uka at the Fermi surface, q„{g,n,m) = 
Uf,- f^{g,n,m) — ^{g,n,m)x^ The rules for the diagrammatic expansion of these ex- 
pectation values in powers of {g^ — 1) were developed by MV, with the result 

00 

d{g, n, m) = / ^(/ - l)^~^Cp(n, m), (4a) 
p=i 

1 /I 2\ 

where n^^ = Ukaid = l,n,m). The functions Cp{n,m) and fp„{k,n,m) can be represented 
by Feynman diagrams corresponding to those of the energy and the Greens function, respec- 
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tively, of a theory. For later convenience we define 

-f e ^ ~ 1^1 /cN 
fla = -riarif,^, foa = Uj^^^ Cq = ^ . (5) 

The diagrams for Cp{n,m) can be obtained from those for fp(j{k,n,m) by connecting the 
two external vertices (see MV): 

Cp(n,m) = --^ ^ra^^/p^(fc,ra,m), p > 1. (6) 
fc 

This equation yields sum rule^Hll for the density of particles inside, n^, or outside of the 
Fermi surface, n^, namely 

n< = n^-n> = ^Y^ nl^Uka, (7) 
fc 

'^f='^r = T^( i d{g,n,m)\. (8) 

It is sometimes useful to remove the diagrams Cp from /pg. and thus define 

hp„{k,n,m) = fpa{k,n,m) - Cp-i{n,m), p>l. (9) 
Note that our definitions in Eqs. (|^) and @ differ slightly from MV. 



B. Canonical transformations 

For the relations to be discussed next the hopping amplitude tij is assumed to be non- 
zero only for hopping between sites i and j on different sublattices A and B. In the next 
chapter we will see, however, that in dimension D = 1 this requirement can be dropped. The 
simplest canonical transformation is the interchange of spin indices t and | which implies 

d{g,n,m) = d{g,n, —m), (10a) 
Ukaig, n, m) = nua{g, n, -m). (10b) 

Furthermore, a particle-hole transformation for both spins yields (MV) 

d{g^ n, m) = d{g, 2 — n, —m) + n — 1, (Ha) 
Ukaig, n,m) = 1- UQ^kaig, 2-n, -m). (lib) 

Here Q is a vector in the first Brillouin zone with e^^^ = ±1 for a lattice vector R E A, 
B, respectively; Q = (vr/a, vr/a, . . . 7r/a) for a hypercubic lattice with spacinga. From 
Eqs. (|10b|) and (|llb|) it follows that the discontinuity at the Fermi surface obeysc2l 

q^{g, n, m) = q^{g, n, -m) = q^{g, 2 - n, m). (12) 

Note that in particular = q^ for n = 1. 

For densities < tt, < 2 the magnetization is in the range \m\ < min(r;,, 2 — n). In view 
of Eqs. (p!OD, (|ll]), and (^2[) we limit ourselves from now ontoO<m<n<l. Therefore a 
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= I will be referred to as "majority spin" and a = J, as "minority spin." Note that the case 
n = m (the fully polarized state without doubly occupied sites) can be obtained from the 
un correlated case (7 = 1. 

Performing a particle-hole transformation for f-operators only,l 

c:^ = (-irc+, c:^ = c,^, as) 

one may derive the following identities for d{g,n,m) and nk^{g,n,m): 

d{g, n, m) = — d{g~^, 1 — m,l — n). (14a) 

nk^{g,n,m) = 1 - nQ^k^{g'^, 1 - m, 1 - n), (14b) 
nkiig,n,m) = Ukiig'^, 1 - m, 1 - n). (14c) 

For the uncorrelated case {g = 1) we have in particular 

nl^in^m) = 1 -n°Q_^^{l-m,l-n), (15a) 
n^|(n, m) = n^i{l — m, 1 — n). (15b) 

The relations in Eq. (|l^) express a property of the starting wave function [Eq. (^)] and can 
also be derived directly from the fact that hopping occurs only between A and B sublattices. 

C. Diagrammatic relations 

We now derive diagrammatic relations for Cp and fpa from the identities in Eqs. (0) and 
(p!5|). The following equations are valid for all p > and < m < n < 1 unless noted 



otherwise. 

1. Double occupation 
From Eq. (K^) we obtain 



/ 1 \ p+i 

dig-\ l-m,l-n) = J^id' - l)^(-l)" - c,+i(l - m, 1 - n). (16) 

p=0 ^ 

Our goal is to equate coefficients of powers of {g"^ ~ !)• To this end we make use of the 
expansion 



^ / r=0 ^ 



We then obtain from Eq. ( p.4a| ) 



Ag. n, m) = - ^(/ - 1)^ (-1)^ Y.\r) ""'^^^^^ -m^-n), (18) 

p=0 r=0 ^ ^ 
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while from Eq. (^a]) we have 



d{g, n, m) = Ci(n, m) + ^(5-^ - 1)^ (cp(n, m) + Cp+i(n, m)). (19) 



(20) 



We are thus led to the relation 

Cp(n,m) +Cp+i(n,m) = (-1)^+^^ P j Cr+i(l - m, 1 - n). 

Now we employ the binomial inversion formul(^ 

ap = {-If J2 Q ^1 all p > 

^ bp= {-ly J2 (^) s all P ^ 0' (21) 

q=0 

which is valid for arbitrary ap and bp. When applied to Eq. ( PD| ) it yields 

Cp(n,m) = i-irYl (l)cril-m,l-n). (22) 

r=0 ^ ^ 

We stress that the relations in Eqs. (|20|) and ( ^2]) are valid in arbitrary dimensions for lattices 
with hopping between A and B sublattices only. In the next section this equation will be 
used to calculate Cp in D = 1. 

2. Momentum distribution 

An analogous procedure is used to derive relations for fp^. We define the abbreviations 

Fp„{k, n, m) = fp+2a{k, n, m) + nl^fp+i„{k, n, m), (23a) 
Fp^{k, n, m) = Fp^{Q — fc, 1 — m, 1 — n), (23b) 
Fp|(fc, n, m) = Fpi{k, 1 — m, 1 — ra), (23c) 

and rewrite Eq. (|4bD as 



^ 00 

Ukaig, n, m) = <^(n, m) + ^(/ - 1)^+^ ^^^^^^ (34) 



~'~ ' p=0 



This expression appears on the left-hand sides of Eqs. ( |14b| ) and ( |14c| ) , while their right-hand 
sides take the form 

1 - nQ-k]ig'\ 1 - m, 1 - n) = n^^(n, m) + ^ Yl (_^2)p+i ^ptC^' (25a) 
^fci(^"\l-"^, 1-^) =<i("',"^) - ■jY^^l^-(Z^iyTr^pi(^'"''"^)- (25b) 
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Expanding {g"^) ^ ^ in powers of {g'^ — 1) [Eq. (0)], comparing coefficients, and combining 
both cases, we find 

Fp^{k,n,m) = -sgn(a)(-l)P^ (^] Fr,{k,n,m). (26) 

r=0 ^ ^ 

For clarity we will from now on label fpa, hp^, and rika with the subscripts < and >, 
depending on whether the momentum lies inside or outside of the Fermi surface. We ffist 
simplify the equations for /p|. Using Eq. (^) together with (|2^) we obtain 

ViT(fc, H = (-1)'' Yl (^) hrMQ - fc, 1 - m, 1 - n), (27) 
valid for all k, relating h^-^ and h^-^. For the minority spin a similar calculation yields 

/< (fc, n, m) = (-1)^ [^"j /< (fc, 1 - m, 1 - n), (28a) 

/>+2i(fc,n,m) = (-l)^X^ (l)f^+2iik, 1 - m, 1 - n), (28b) 

i. e., there is a relation between the /p| for momenta inside of the Fermi surface, and another 
relation for momenta outside of the Fermi surface, each linking the cases n + m < 1 and 
n + m > 1. Note that the relations in Eqs. (P7|)-(PB|) are valid in arbitrary dimensions for 
lattices with hopping between A and B sublattices only. 



III. ANALYTIC EVALUATION IN D = 1 

In the remainder of the paper we consider the GWF for the Hubbard chain with a 
symmetric dispersion = e_fc that increases monotonically with \k\. This implies that 
nka = n^ka] therefore we only consider k > 0. For convenience we assume that the zero of 
energy is chosen such that the center of the band, e = 2 dkek, is zero. 

For such dispersions the free Fermi sea described by the starting wave function |$o) is 
centered around k = and is simply connected, i. e., = Q{k-Fa — \k\) is a step function. 
(We follow the convention of MV of measuring k in units of 27r/a, where a is the lattice 
spacing; the first Brillouin zone is the interval [—1/2; 1/2], the reciprocal lattice vectors K 
are integers, and the nesting vector Q = 1/2.) For this Fermi surface topology the Fermi 
momentum, kp^ = n„/2 = {n + sgn((j)m)/4, only depends on the particle density i. e., 
the particular form of the dispersion is irrelevant. In general this simplification occurs only 
in dimension D = 1; higher-dimensional tight-binding dispersions usually do not depend 
only on although this symmetry can be artificially imposed to allow the construction of 
a dispersion from a given density of states.ii 

Since the values of the diagrams Cp and fpa are completely determined by the function 
n^^ (which is independent of the dispersion as long as et is increasing with |A;|) the relations 
derived in Sec. [II B| , and hence the relations in Eqs. (pOD, (|2^) , (pT]), (|^), are valid for all 
increasing and symmetric dispersions m D = 1. 



7 



The following analytic calculation of GWF expectation values with magnetization m 7^ 
is based on the corresponding calculation for m = by MV. The calculation for m = 
was made possible by exploiting the relations following from canonical transformations, the 
polynomial form of the diagrams and their continuity as functions of k and n, and an analysis 
of the contribution of Umklapp processes. We will now use very similar methods to express 
the double occupation d{g,n,m) and the momentum distribution nk„{g,n,m) in terms of 
the known quantities for m = 0. 

In Section A we review the results of MV and present new closed formulas for n^o- for zero 
magnetization. For the magnetic case the double occupation, the momentum distribution, 
and its discontinuity at the Fermi surface are calculated in Sections B, C, and D. The 
variational energy is evaluated in Section E. 



A. Zero magnetization 

For the non-magnetic case the diagrams Cp and hp^j were already calculated by MV. For 
m = 0, n<l, 0<k< 1/2, the results may be summarized as 

Cp(n,0) = 'jpU^^'^, 



fpa{k,n,0) = hp„{k,n,0) + Cp_i(ra, 0) 



ioi k <j 

for f < A; < min(f ) 
for f < A; < i 
for 1-f <k<\ 



(29) 
(30) 

p>2. (31) 



Here and below 7^ = (— l)^+^/[2(p + 1)], and Rp{k) and Qp{k) are certain polynomials of 
degree p in k. Note that at k-p = n/4 both fp„ and hp^j are discontinuous, and hp„ = for 
k > 3n/A. The contribution with momentum 1 — in Eq. (0) is due to Umklapp processes. 
From the diagrammatic series in Eq. (S) MV's result for the double occupation is 



d{g,n, 0) 



9' 



2(1-9')' 

and the momentum distribution n^o- is given by 



(-ln[l-(l-5>]-(l-^7». 



1- 9 n 



n,A9,n,0) = l~-^^ 2 

1 _ n 

nkai9,n,0) = 



+ 



9' 



;i + 
1 



[7^o( 



ik 



9»-ll 



[-iln[l-(l-^2 



1 — g n 
' TT^ 2 + {l+gy 

+ Qo{^,il-9')n) + QoC-^^^, 



n 



9 )n 



for k < k-p and k > kp, respectively. Here we introduced the generating functions 



(32) 



(33) 



(34) 



TZo{x,z) = ^R, 



p=0 

00 



p=0 



(35) 
(36) 
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with the convention that TZq{x,z) and Qo{x,z) are zero for \x\ > 1. 

In MV the coefficients of the polynomials Rp{k) and Qp{k) had to be calculated recursively 
to obtain the momentum distribution, and convergence was problematic for (1 — g'^)n close 
to 1. In Appendix ^ we show how the recursion equations for Rp{k) can in fact be solved 
in closed form [Eqs. ( |A5| )-([AE|)], while Qp{k) can be expressed in terms of them and their 
integrals. The generating functions are calculated as 



v/(2 - zf - [xzf \^{2-zY -[xzfr ^ ^ 

Qo(a;, z) = Wo{x, z) + ^ [(1 - x) 7^o(a;, z) + 7^l(a;, z)], (38) 

for |x| < 1. Here K{k) = /^^^^(l — k"^ sin^ (j))^^^^d(j) is the complete elliptic integral of the 
ffist kind, yVo{x,z) is an auxiliary function, 

yVo{x,z) = ^^no{x,z) + ^^ni{x,z) + ^^^-^ni{x,z), (39) 

where the dot indicates a partial derivative with respect to second argument, and 7lj{x,z) 
is the repeated integral of 7lo{x, z), defined by (j > 0) 

TZj+iix, z) = dx' TZjix', z) = ^ J^ dx' {x - x'Y 7^o(x , z). (40) 

Below we will also need the following integral, 



/X 
dx' Qoix',z) 



(1 — z)(x — 1) ^ , . 5z — 4 ^ , , z(z — l) d ^ , , 

= ^ 7^l(x, z) + 7^2(x, z) + ^ ^ ^ — 7^2(x, z). (4i) 

In Appendix ^ we provide an explicit expression for 7lj{x,z) [Eqs. ( [A10D -( |ATT| )], as well 
as other relations. Here we note in particular the following functional relations, which are 
obtained from Eqs. (|37D and (^), 

7^,■(x,^) = (l-2;)7^,(x,2;), (42) 
Q^{x,j^,) = W,{x,z). (43) 

These expressions analytically continue 7lj{x,z) and Qj{x,z) to \z\ > 1. 



B. Non-zero magnetization: double occupation 

To calculate the double occupation one needs the diagrams Cp. Using the methods of 
MV one can show that for m > Umklapp processes occur in the graphs for Cp{n,m) 
only if 72 > 1. This implies that for n < 1 Cp{n,m) is a homogeneous function of kp^ 
and kpi of order p + 1, since every contributing graph contains p + 1 free momenta. Since 
kpa = {n + sgn((T)m)/4 it follows that Cp{n,m) can be written as a linear combination of 
terms nln^^^ *" with < r < p + 1, i. e., Cp{n,m)/n^^^ is a polynomial in m/n of degree 



9 



p+1. This is also the case for n = 1 since d{g, n, m) is continuous at n = 1 [see Eq. ( p.la| )]. 
Therefore we can write the polynomial simply as Cp{n,m)/n^~^^ = Cp{l,m/n) for n < 1. For 
n = 1, however, we can calculate Cp(l,m) from Eq. ( p2D and Eq. (pQl): 



Cp{l,m) = {-iyJ2 -m,0) 

r=0 ^ ^ 

= 7p(l-mP+'). 
The polynomial structure of Cp{n,m) then implies 



(44) 



Cp(n,m)=nP+ic,(l,f) = 7p( 



p+i 



= Cp{n,0) - Cp(m, 0). 
Summation of the series in Eq. ( ^a]) yields the simple result 

d{g, n, m) = d{g, n, 0) — d{g, m, 0). 



(45) 



(46) 



It is remarkable that the double occupation at density n and magnetization m is obtained 
as the difference between the double occupation without magnetization at density n and 
density m. With the closed form of d{g, n, 0) taken from Eq. ( p2D we finally obtain for the 
double occupation, valid for < m < < 1, 



d{g,n,m) 



9' 



2(1 -g 



2\2 



1 l-(l-^^)m ^ 

^^1 n ^\ ^ ~^ ){n-m) 

1 - {^~9 )n 



(47) 



The double occupation is shown in Fig. |T| for various parameter values. In the limit of strong 
correlation [g —>■ 0) it behaves as 



,2 r 



d{g, n < l,m < n) = — 
d{g,n = l,m < 1) = 



1 - m ' 
m [n — m) 



n 



0(9'), 



■ 1 ln(l - m) - (1 
In - H ^ ^ 
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+ 0{gHng), 



(48a) 
(48b) 



i. e., the double occupation is non-analytic in the limit n ^ 1, g ^ 0, asm the non-magnetic 
case (see MV). 



C. Non-zero magnetization: momentum distribution 

To calculate the momentum distribution n^a one needs the diagrams /pg- or hpo-. In the 
following subsections we distinguish the cases of majority and minority spin, as well as 
whether k lies inside or outside of the Fermi surface. 

For finite magnetization, m > 0, we make the following observations, similar in spirit to 
those of MV for the non- magnetic case. As a function of k, fpa{k, n, m) and hp(j{k, n, m) are 
discontinuous at = /cpo-, since the one-particle irreducible graphs contain a factor n°^. For 
k > /cpcr + momentum conservation at the outer vertices of hp„ cannot be fulfilled, so 
that in this case hpcr{k, n, m) = 0. For k > 1 — k-pa — 2kp^ (> kpcr) Umklapp processes occur 
and yield an additional contribution to hp^ of normal processes with external momentum 
1 — k. In the absence of Umklapp processes fpa and hpa are homogeneous functions of /cpt 
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and kpi of order p, since every contributing graph contains p free momenta. Then fpa/n^ 
and hpa/n^ are polynomials in k/n and m/n of degree p, and due to momentum conservation 
at outer vertices different polynomials occur depending on whether k is larger or smaller 
than -^{2k-pa — k-pa)- Furthermore, different polynomials for h^Jn^ occur also depending on 
how k compares to 3/cf| (see subsection 4 below and Appendix |^). 



i. Majority spins inside of the Fermi surface: n^^{g,n,m) 



For < k < kp^ = (ra + m)/4 and n = 1 the diagrams h^.^ can be obtained from Eq. (27) 
in terms of the known functions h^-^ of the non-magnetic case [Eq. (|3TD ]: 



hp+i^{k, l,m) 




a-k,l-m, 0), 



(49) 



for <k<^ 
for < A; < 



3m— 1 



(50) 



1— 3m 



Since Umklapp processes do not occur for k < kp^, hp^{k,n,m)/n^ is a polynomial in k/n 
and m/n. Hence we can obtain h^^^^/n^^^ by replacing khj k/n and m by m/n in Eq. (pOD. 
This yields 



h<^,,{k,n,m) = i-irj2(l) 



{n — my n 



Q 

X { 

[Q 



r+V 



2 t\ 

n—m ' 



for <k<^ 



for < A; < 
+ Qr+i(S)] forO<A;< 



3m— n 



ji- 3m 



(51) 



Note that hp^^.^{k,n,m) is continuous for all k < kp^- due to (5p(3/4) = (see MV). The 
normal processes contribute differently depending on how k compares to ±{2kp^ — kp^), as 
expected. 

Since h^-^ is available for arbitrary orders of p, the series in Eq. (^) can now be summed, 
with the result, valid for 0<m<r2<l, 0<A;<(n + m)/4, 



nk^{9,n,m) 



[I + 9)' 



n — m 



- d{g,n,m) 



VA) ; 2m-4k (l-g^)(n-m) 
•^^0 1 n-m ' l-(l-g2)m 



Wo 



2m+4fc (l-g^)(ra-m) 
n—m ' 1 — (1— g2)m 



(52) 



where the functional relation in Eq. ( ^31) was used. For m = this result reduces to Eq. (^ 
by virtue of Eq. (^). 
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2. Majority spins outside of the Fermi surface: n^-^{g,n,m) 

For (n + m)/4< k < 1/2 we can deduce h^-^ at once by inserting h^^ into Eq. 



r=l ^ ^ 



{n-my+\-my 



fe— — 



X { 



for m < A; < min 
for 3m < A; < 1 



' 3n—m 
, 4 5 



-1 1. m 

' ^ 2 M -f^-v 1 _ 3n— m 



-)] for 1- 



4 <k<l 



3n—m ^ 



(53) 



which is continuous for all k > kpp In this sector Umklapp processes occur for k > 1 — 
2kpi — kp^ and contribute to /iS like normal processes with external momentum 1 — k. 

Summing the series in Eq. or applying Eq. ( |14b|) to Eq. (|52|) , yields the result, valid 
for < m < n < 1, {n + m)/A < k < 1/2, 



d{g,n,m) 



4fc-2n (l-g^)(n-m) 
n—m ' l—{l—g'^)m 



4(l-A:)-2n (l-g^)(ft-m) 
n—m ' 1 — (1— g^)m 



(54) 



which alternatively can be derived directly from Eqs. ( |14b|) , (|5^), and (0). Note also that 
for m = this result reduces to Eq. (0). The momentum distribution nki{g, n, m) has thus 
been determined in the entire parameter range. It is shown for g = 0.1 and densities n = 1 
and 0.8 in Fig. 0a and Ha. 



3. Minority spins inside of the Fermi surface: n'j:^{g,n,m) 

For < k < kpi = {n — m) /A and n = 1 the diagrams f^^ are given in terms of the known 
functions of the non-magnetic case [Eq. (|31|)] according to Eq. ( p8a| ): 

f<ik, 1, m) = i-iy X: h (1 - mYRrij^). (55) 

T- = ^ ^ 

Due to the absence of Umklapp processes for k < kpi, fp^{k, n, m)/n^ is a polynomial in k/n 
and m/n. Therefore fpjn^ is given by Eq. (|55D with k replaced by k/n and m replaced by 
m/n. We then have 




where the second equation was obtained by inserting Eq. ([56| ) in Eq. ( |28a|) ; it is equivalent 
to the functional relation in Eq. (^2]). Note that for m = Eq. ( ^Tf ) reduces to (pl|). 
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We use Eq. ([57|) in Eq. (|4bD to find the following expression for n^^, valid for < m < 

n < 1, < k < (n - m) / A, 



n^^{g,n,m) = 1 



1 — g n + m 



+ 



1 + g 2 (1 + ^) 



7^, 



4fc (l-g^)(n-m) 
1 n—m ' 1— (1— g'^)??! 

1 — (1 — (yf^)m 



- 1 



(58) 



This expression reduces to Eq. (^3|) for m = 0. 



^. Minority spins outside of the Fermi surface: n-^^{g,n,m) 



Finally we consider the case (n — m)/4 < k < 1/2, for which the calculation of n-^^ is 
somewhat more complicated. We begin with the special case n = 1, for which the momen- 
tum distribution can be determined immediately from Eqs. ( p.4c[ ) , ( |10b| ), and ( |14b[ ), which 
combine to give 



^kiiaA^m) = n^j^ig \l - m,0) = n^^ig \l - m,0) = 1 - nf_^^{g,l,m) 



(59) 



Inserting Eq. (p2D we arrive at 



+ 



9' 



- d{g, l,m) 



f 2m+ik-2 (l-g^)(l-m) 
I 1— m ' 1— (1— g^)m 



I / 2m-ik+2 (l-g^)(l-m) 



(60) 



Next we consider arbitrary density and magnetization, < m < n < 1. We make use of the 
following relation, which follows from Eqs. ( P8b|) and (pO|), 



h>^2iik,n,m) = -2cp+i{n,m) - Cp{n,m) + (-1)^^ P j 



(^; 1 — "^j 1 — ""-)• (61) 



The diagrams h^^ that appear in this equation can be written in terms of the contribution 
of normal processes, Np, as 



' k m,\ 



hpi{k,n,m) 



for ^^^j^ < k < min 
for 3H±!ii < A; < I 



' 3n+m 
, 4 5 



3n+m 



(62) 



n^'[iVp(^,^)+iV,(i^,^)] forl-i^<A;<l 



In Appendix ^ we show that depending on k and m, the function Np{k, m) is given piecewise 
by four polynomials in k and m of order p. The explicit determination of these polynomials 
is quite involved; it is presented in Appendix 0. The final result for n^^, valid for < m < 

n <1, [n — m)/A < k < 1/2, can be written as 



n^iig, n, m) = f' d{g, n, m) + 



(1 + , 



[N{g, k,n,m) + N{g,l — k,n,m)], (63) 
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where N{g, k, n, m) is given by 



N{g, k, n, m) 



' N^'^\g,k,n,m) 
N'^'^\g, k, n, m) 
N'^^\g, k, n, m) 
N^^\g, k, n, m) 




for M < A; < min(^^, ^^^) 
for m < f and ^ <k< 
for m > f and <k< 
for max( 



^) < A; < ^ 



for 



3n+m 



< k 



with 

N'^^\g,k,n,m) 



(1 - (1 - g^)m) (3n - 3m - 4A;) 



+ 



+ 



2{n — m) 
[n — 3m) (1 — g"^ 



Qo 



4fc-2n+2m (l-g^)(n-m) 
n—m ' 1 — (1— g*^)rrt 



4 



^ Q / 4fc-2n+2m (l-g^)(n-m) 
^ \ n—m ' 1— (1— (j|2)m 



[n 



m) {I - g^) {I - {I - g^)n) ■ Uk-2n+2m {l-g^)(n~m) 

2(1- (I- g'^)m) V " ' 



iV(^)(c/,fc,ri,m) 



+ N^^\g,k,n,m), 
1 - (1 - 5'^)m) [Ak-n - 3m) 



2 (ra — m) 
(n - 3m) (1 - ^2) + 2 



Qo 



4fc-2n-2m (l-g^)(n-m) 
n—m ' 1 — (1— g^)m 



4fc-2ra-2m (l-g"')(n~m) 
n—m 



N'^^\g,k,n,m) 
N^'^\g,k,n,m) 



(n — m)(l — (yf^) 



^ y n—m ' 1— (1— g-^)™ 

m) (1 - g^) (1 - (1 - g^)n) l^ 4k-2n-2m (l-g^)(n-m) 

1 — g'^)m 



2{l-{l-g^)m) 



l-(l-/_)^j^l 



2 2 1 - (1 - g^)n 

N'-^^ [g, k, n, m) - N'-^^ [g, k, n, m) + A^^*^) [g, k, n, m 



n—m ' 1— (1— g^)m 
2^ 

)iib 



(64) 



(65a) 



(65b) 

(65c) 
(65d) 



and the dot again denotes derivative with respect to second argument. It can be checked that 
Eq. (|63| ) indeed reduces to (|60| ) for m = 0. Thus the momentum distribution nki{g,n,m) 
has been determined for all parameters. It is shown for g = 0.1 and densities n = 1 and 0.8 
in Fig. |b and |b. 

The calculation of the correlated momentum distributions n^o- is now complete. We 
remark that they are continuous functions of k, except at kp„, and are also continuous in n 
and m for fixed k. It can be checked that they obey the sum rule in Eq. (H). 



D. Discontinuity of the momentum distribution at the Fermi surface 

It suffices to calculate the discontinuity of Uka at the Fermi surface for < m < n < 1 
[see Eq. (^2])]. From our previous results we obtain the left and right limit of the momentum 
distribution at the Fermi vector as 

n^^Ja, n, m) = n, + (66b) 
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where the abbreviation is defined as 



= y/[l-(l-(^2)^][l_(l_^2)^]-sgnM. (67) 

For the discontinuity at the a-spin Fermi surface we thus obtain 

g.(^,-,m)=(^^^. (68) 

It follows that vanishes only for a half-filled band without double occupation (rz = 1 and 
g = 0); in this case there is exactly one particle at each site so that = 1/2 for all k. 
Note also that = if n = 1 (or, trivially, if m = or g = 1). We plot q^^ for n = 1 and 
n = 0.8 in Fig. |. 



E. Energy expectation value 

For any symmetric dispersion e^, monotonically increasing with \k\, we can now calcu- 
late the energy expectation value per site, Eq = {H)q/L, of the one-dimensional Hubbard 
Hamiltionian (|I]), which is then minimized w.r.t. g to find the optimal variational energy, 

EQ{n,m,U) = min EQ{g,n,m,U), (69) 

0<5<1 
1/2 

EG{g,n,m,U) = 2 dk ek ^nka{g,n,m) + U d{g,n,m). (70) 



Note that it follows from Eq. ( |llb|) that the total kinetic energy for dispersion at density 
n > 1 can be calculated from the dispersion —ei/2-k at density 2 — n using the formulas for 
Uka and d that are valid below half-filling. 

For the Hubbard chain with nearest-neighbor hopping t the dispersion relation in our 
notation is ek = —2tcos{27ik). We assume t > without loss of generality, so that the 
dispersion is increasing with \k\ and our results for the Gutzwiller expectation values apply. 
The optimal variational energy for this system is shown in Fig. |^ for densities n = 1 and 
0.8 for various magnetizations. Note that at half-filling no Brinkman-Rice metal-insulator 
transition occurs at any finite U; i. e., g = is the optimal variational parameter only for U 
= oo. The variational result for the ground-state magnetization is determined in the next 
section. 



IV. MAGNETIC PHASE DIAGRAM OF THE HUBBARD CHAIN 

In this section we determine the instability towards ferromagnetism for the Hubbard 
chain with nearest-neighbor hopping. Currently only homogeneous paramagnetic and fer- 
romagnetic phases can be investigated analytically with the GWF in D = 1; hence we do 
not consider antiferromagnetism or other broken symmetries. We begin by examining the 
energy for the special cases of the paramagnetic state (i. e., zero magnetization) and the 



15 



fully polarized state. The latter contains the minimum number of doubly occupied sites and 
is an eigenstate of H ^ with eigenvalue 



i?Fp(^, U) = Eq{1, n, min(n, 2 — n),U) 

eo{n) for < n < 1 

eo(n- 1) + 1) for 1 < n < 2 ' ^ ^ 

where eo{n„) = 2 ^^''^'^dkek is the kinetic energy of one spin species for the uncorrelated 
state. For the case of nearest-neighbor hopping we have eo(no-) = — 2t sin(7rno-)/7r, t > 0. 

In Fig^ the exact ground-state energy, E{n, m = 0,U), obtained from the Bethe ansatz 
solution,BJ is compared to the Gutzwiller energy for zero and maximal polarization, at various 
densities. We also show the energy of the variational Hartree-Fock theory, E}ip{n, m = 0, f/); 
it is contained as a special case in the results for the GWF, 

EiiF{n, m, U) = Eg{1, n, m, U) 

n + m\ fn — m\ U , ^ 2\ 



As expected, the (^-optimized GWF significantly improves upon Hartree-Fock theory but 
overestimates E at large uM Since all spin configurations are degenerate for U = oo, the 
exact ground-state energy coincides with Epp in this case; therefore the Gutzwiller energy 
Eq necessarily crosses the value E-pp at some finite value of the interaction Uc (except for 
n = 1). The existence of a finite critical interaction Uc above which the GWF predicts a 
ferromagnetic ground state is in contrast to the Lieb-Mattis theorem,!! which states that 
m = for the exact ground state (i. e., Uc = oo). The reason for this overestimation of 
the instability of the paramagnetic state lies in the simple structure of the GWF, which 
controls only local correlations and cannot describe the special correlated behavior in D = 1 
microscopically. 

The preceeding discussion only compared the variational energies for zero and full polar- 
ization. From our new results for the ferromagnetic GWF we can also study the stability 
of partially polarized ferromagnetic states. We first consider Hartree-Fock theory. A simple 
calculation shows that it predicts a fully polarized ground state for U > U^^{n), where 

rrHF. ^ _ /l6sin(7™/2) [1 - cos(7rn/2)]/(7rn2) for < n < 1 

^''^"\f/H^(2-n) forl<n<2- ^^^^ 

This critical interaction U^^ is smaller than that derived from the Stoner criterion, U^^{n) = 
l/N{en/i) = 27r sin(7rn/2), where N{e) is the density of states. Note that Ebf as a function 
of m never develops a local minimum at m 7^ 0. On the other hand, a maximum at finite 
m occurs for U > 2(1 — cos(7rn))/ min(n, 2 — n), which leads to a global minimum at full 
polarization already for U > U^^. The Stoner criterion, which merely signals a negative 
curvature of E^p at m = and does not take into account a finite magnetization, is thus 
irrelevant for the Hubbard chain with nearest-neighbor hopping. 

For the Gutzwiller wave function we find that Eq as a function of m at fixed U devel- 
opments local extrema and global minima in a qualitatively similar fashion to -Erf- As a 
consequence the GWF also describes a discontinuous transition from the paramagnetic state 
to a state with full polarization at U = Uc{n). This critical interaction Uc is shown in Fig. |^. 
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Compared to Hartree-Fock theory we find agreement in the limit of small n. However, at 
intermediate densities the GWF predicts a significantly reduced ferromagnetic region. In 
particular for n — >• 1 we have Uc ^ oo, as expected from the previous discussion (see also 
Fig. ||a). Thus, in contrast to Hartree-Fock theory, the GWF does not exhibit a spurious 
ferromagnetic transition at half-filling, since it is able to avoid double occupation not only 
through a ferromagnetic polarization, but also by decreasing the variational parameter g. 
Away from half-filling, however, the GWF predicts ferromagnetism for sufficiently large U, 
in contrast to the exact solution for the Hubbard chain. 

V. CONCLUSION 

In this paper we presented general diagrammatic relations for the expectation values 
of the D-dimensional Hubbard model in terms of the Gutzwiller wave function (GWF) at 
non-zero magnetization m. In D = 1 explicit, approximation-free evaluations of the double 
occupation d{g,n,m) and the momentum distribution nkcr{g,n,m) were made possible by 
exploiting (i) relations for the Feynman diagrams for d and Uka, derived from canonical spin 
and particle-hole transformations, (ii) the polynomial form of the diagrams in powers of 
k, n, and m, and (iii) an analysis of the contributions of normal and Umklapp processes. 
In this way the calculation of d{g,n,m) and nk„{g,n,m) was reduced to that for m = 
0. Furthermore, new closed expressions for the momentum distribution n^a were derived, 
facilitating numerical evaluation. 

The functions d{g,n,m) and nk(T{g,n,m) in D = 1 are qualitatively similar to those for 
m = 0. The discontinuity q^ig, n, m) of the momentum distribution at the Fermi energy was 
also calculated explicitly. It is always finite, except for the half-filled band without double 
occupation {g = 0, n = 1, m = 0) in which case the electrons are trivially localized. In all 
other cases the GWF describes a ferromagnetic Fermi liquid. 

Analysis of the Gutzwiller variational energy for the Hubbard chain with nearest-neighbor 
hopping shows that the GWF predicts a fully polarized ferromagnetic state at large enough 
U and away from half-filling, in contrast to the Lieb-Mattis theorem.il This exemplifies 
once more the peculiarities of the GWF which controls correlations between the electrons 
only globally through the local Hubbard interaction. While the GWF is an excellent wave 
function for the one-dimensional Heisenberg model (at least for m = 0)p since this only 
involves spin correlations between localized spins, it is not a very good wave function for 
the one-dimensional Hubbard model at large U and n ^ 1 since it does not describe density 
correlations well in this case.a As a consequence the GWF cannot reproduce all character- 
istics of the one-dimensional system. This is also apparent from the finite discontinuity of 
the momentum distribution at the Fermi surface, which is, in fact, continuous for Luttinger 
liquids such as the one-dimensional Hubbard model. On the other hand the ferromagnetic 
GWF represents a trial state for partially polarized, itinerant electrons and may thus be 
regarded as an effective, non-perturbative description of a ferromagnetic Fermi liquid. 

In view of the considerable technical complications involved in the present calculations 
it is not clear whether it will be possible to compute correlation functions with the GWF 
for m 7^ 0. Since the calculation of the spin-spin correlation functioni for m = helped to 
gain considerable insight into the properties of Heisenb erg-type modelsj^'i a corresponding 
result for m 7^ would be helpful for a better understanding of one-dimensional Heisenberg 
models in a magnetic field. 
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APPENDIX A: CALCULATION OF POLYNOMIALS FOR THE NON-MAG- 
NETIC CASE 



In this appendix we describe the derivation of closed expressions for the polynomials Rp{k) 
and Qp{k) that appear in Sec. |ITT| . By eliminating Qp{k) from MV's recursion formulas we 
obtain 



p2 Rp{k) + p{p + 1) = k {{2p - 1) R'pik) + 2p 

R;{k) + R;^,{k) 



Furthermore, the polynomials Qp{k) can be expressed in terms of Rp{k) as 

Q'p+iif" + 1) = + 1) Rpik) - 2(P + 1) Rp+i{k) 

+ (2A;-i) {R;{k) + R'^^,{k)), p>0, 

together with Qp(|) = 0. We define Ro{k) = 1, Qo(^) = 0. 

A closed form for Rp{k) is obtained as follows. Using Eqs. (|5 
Eq. Q to 



p>0. 



(Al) 



(A2) 

we can reduce 



p{p + 1) Rp{k) + 2k R'{k) + {e - ^) R''{k) = p" Rp-i{k), p>0. 



(A3) 



This is essentially the differential equation of the Legendre polynomials except for the in- 
homogeneity on the right-hand side. From an expansion in Legendre polynomials Pn{x) we 
thus obtain 



L|pJ 

Rpik) = J2 



j=0 



{-iyp\\4j + i) 
(p - 2jy.{p + 2j + iy\ J 



IN 2 
2 



P2,(4fc), p>0, 



(A4) 



which after some calculation yields Eq. (pTl). From TZq{x,z) we then obtain the explicit 
expressions in terms of a (terminating) hypergeometric function. 



j=0 



4^+P{p ~ 2j)\{p - 



p 



5i^2 + i I - P. 1; 1 - 



p 



2 \ ^^F,{-p,-p,\-\ 



p 



JJ \ J J \P-J 



1 . 4fc-l N ( ik+l \P 
t'-' Ak-^\l \ 2 ) 

\-^k\ fl^^kV"^ 



(A5) 

(A6) 
(A7) 
(A8) 
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Furthermore an integration by parts of Eq. leads to the expression for Qo{x, z) in terms 
of 7lj{x, z) shown in Eq. (pSl). 

By using a hypergeometric identity to rewrite Eq. ( pT]) as 

' <1, (A9) 



(a;2-l)z2\ 


(l-a;2)22 


4(1-2) ^' 


4(1-2) 



^o(^, 2;) — , 2F1 I 2, 2' 

Vl - ^ 

integrating term-wise w.r.t. x, and again using several hypergeometric identities we obtain 
the following explicit expression for 7lj{x,z), 

UJx, z) = y ^-^-^ ' ^ I (AlO) 

where the Rp\x) are polynomials of degree p, 

= (7) (7+'^.) M-P^hrA + P + 3.\+3M)^^ + ^y- (All) 

The functions TZj{x, z) may be evaluated via the series ( |A10|) for not too large values of z. 
Alternatively, the integration in Eq. (|40|) can be performed numerically. 

Finally, we note a few special values. From MV's polynomial relations we obtain 

Qo{x, z) + Qo{-x, z) = l-{l-z) TZoix, z) + ^ ln(l - z). (A12) 

Together with Eq. (^)-(|39|) we find in particular 

7^o(±l,^) = ^=^, Qo(l,^)=0, (A13) 
\/l — z 



Qo{~l,z) = l-VT^+-lnil-z). (AM) 
Furthermore the sum rule for particles outside of the Fermi surface [Eq. (|])] implies 

ni{-l,z) = 27^l(0,z) = - ln(l -2), (A15) 

z 

Qi(-1, z) = -l + ^ ln(l - z). (A16) 

APPENDIX B: POLYNOMIAL STRUCTURE OF MINORITY SPIN DIA- 
GRAMS OUTSIDE OF THE FERMI SURFACE 

In this appendix we derive the polynomial structure of the diagrams /i^, for which a 



lot of cases must be distinguished. As discussed at the beginning of Sec. [Ill Q , different 
polynomials occur in the contribution of normal processes depending on the relation of k 
to 2/c| — ki and 2k^ + k^. Furthermore, due to a certain class of diagrams, shown in Fig. |^, 
different polynomials for h^Jn^ may in principle occur also at Sfepi, S/cf;, and all higher 
odd multiples of /cpj,; however some simplification takes place, as we will show below. 
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Case ^4: < m < n/3. Here the momenta 2k^ 
follows, 



ki, 2k^ + k\ , 3fc|, etc., are ordered as 



n + 3m 3(n — m) 3n + m 

< — — : < : < 



5{n 



m] 



4-4 4 - 4 ' 

so that the contribution of normal processes to h^^, as it appears in Eq. 
in terms of three polynomials Ap\k,m) of order p: 



(Bl) 

can be written 



Np{k,0 < m < |; 



'Ai^\k,m) 


for 


l-m ^ u ^ l+3m 
4 - - 4 


A'^\k,m) 


for 


l+3m ^ U ^ 3(l-m) 
4 ^ - 4 


' A'^\k,m) 


for 


3(l-m) ^ u ^ 3+m 
4 - - 4 


.0 


for 


3+m < ^ 



(B2) 



\i k > \ — (3n + m)/4 Umklapp processes contribute to with momentum 1 — fc, and 
several cases must be distinguished to determine the appropriate polynomials. We find for 
n + m < 1 and m > 3n — 2 



/ipi(A;,n,0<m< f) 



/If 
/If 





^ n ' n ■ 
n ' n ^ 



for 



for < A; < 

3{n— m) 



< A; < 

3(n— m) 
4 



for 



for 1 
for 



3n+m < ^ < 1 



3n+m 



(B3) 



while for n + m < 1 and m < 3n — 2 

(2)/fc 



hp^{k, n,0 < m 
nP 



< 



for 



< A; < 



for < A: < 1- 



3n+m 



for 1- 



371+m 



< A; < 



1 I 3(n— m) 



^ ^n^ n ' ^ ^ n ' n ' 

4'^(-,-) + /4?^(— ,-) 



for i_3(^ < i 
for fcl^ < A; < i 



(B4) 



whereas for n + m > 1 



/i>(A;,n,0 < m < f) 





m " 
n ' 






for "-"^ < A; < 1 


^p U' 


m " 
n ' 


+ 4^^(¥' 


m \ 
n ' 


for 1 < A; < "+4*"^ 


m, " 
n ' 




m \ 
n ' 


for "+^3™ <k<\ 1 


^p U' 


m " 
n ' 




m \ 
n ' 


for 1 < A; < 1 


^(3)/fc 
In' 


m ' 
n ' 


+ 4^^(^, 


m \ 
n ' 


for < A; < i 



Case B: n/3 < m < n/2. Here, on the other hand, we have the ordering 
3(?7, — m) ^ n + 3m ^ 5{n — m) ^ 3n + m ^ 7{n — m) 



(B5) 



(B6) 



20 



and the contributions of normal processes to h^^ are now given by four polynomials Bp\k, m) 
of order p: 



Np{k, i < m < i) = < 



B^p\k,m) for ^ ^ ^ ^ 

Bf\k, m) for < A; < 

Bf\k,m) fori^<A;<^M 

5i^)(A;, m) for < A; < ^ 

for ^ < A; 



1— m 



< A; < 



l+3m 



(B7) 



Let us determine the Umklapp process contributions for the region 1 — n<m<n — 1/2, 



hl^{k,n,l<m<l) 



nP 





'k 
^ n ' 


m " 
n ' 


^Bf\ 






for < A; < 1 3'^+'" 




'fc 
^ n ' 


m " 
n ' 


i-fc 

n ' 


m N 
n ' 


for 1 < A; < 1 5(n-m) 




'fc 
^ n ' 


m " 
n ' 


^Bf[ 


1-fc 
n ' 


m N 
n ' 


for 1 '("^ ™) < A; < '("-'"^ 




'fc 
^ n ' 


m " 
n ' 




1-fc 
n ' 


m \ 
n ' 


for '("p) < A; < 1 r+^^m i| 




'k 
^ n ' 


m " 
n / 




1-fc 
n ' 


m \ 
n 1 


for < A; < 1 




'k 
n ' 


m " 
n ' 


+ 5f ( 


1-fc 
n ' 


m \ 
n 1 


for 1 "+43™ < A; < 2 



(B8) 



Now we connect the polynomials with one another via Eq. (0), which performs the trans- 
formation n ^ 1 — m and m \ — n. This provides a link between Eqs. ( P^ ) and 
the right-hand sides of which are thus related, line by line, via (0). Similarly, Eqs. 
and ( [B3| ) are connected, which reveals that the distinction at momentum 1 — 5(?7, — m)/4 
is in fact absent in Eqs. (P5|). Hence we find B^\k,m) = Bp^\k,m), and also, by 

?(3)/ 



m 



comparison with the first transformation, B'^\k^m) = Ap^' {k,m), B^\k^m 
By inspecting the region 2r2 — 1 < m < 1 — raof case B, we find a similar connection to larger 
m (i. e.. Case C: n/2 < m < 3n/5), and use of Eq. (|61]) shows that the new distinction at 
momentum 7{n — m)/4 disappears in a similar fashion, and that all its polynomials likewise 
reduce to the ones above. It is not difficult to see that this simplification repeats for larger 
magnetization, i. e., also for all 3?t,/5 < m < n. We conclude that new polynomials at 
higher odd multiples of [n — m)/A are ruled out by the symmetry of the diagrams; only the 
distinction at momentum 3{n — m)/4 survives. Therefore the contribution ( |62D of normal 
processes to h^, can finally be written as 



A^^\k,m). 



Np{k, m) 



' A^p\k,m) 
Af\k,m) 
Bf\k,m) 
Af\k,m) 




for 



1— m 



< k < min 



l+3m 3(l-m) - 



for m < I and < A; < ^^i^ 

for m > i and ^^i^ <k< 

for max(i^, < A; < ^ 



(B9) 



for 



3+m 



< k 



1. e 

pendix |C[ 



a total of four polynomials are needed to describe n^^] they are determined in Ap- 
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APPENDIX C: CALCULATION OF POLYNOMIALS FOR MINORITY SPINS 
OUTSIDE OF THE FERMI SURFACE 



In this appendix we determine the polynomials that appear on the right-hand side of 
Eq. ( |B9[ ) ■ First let us examine A'^\k,m). For m = we immediately obtain A^\k,0) = 
Qp{k) by comparison with Eq. (|31|). Furthermore we can derive its behavior for small m 
from the equation 

hpi{k,n,m) = h^^{k,n,-m) , (CI) 

which is a simple consequence of Eq. (|10b|) . Similar to MV it can be shown that hp^{k, n, m), 
when regarded as a function of n, has two continuous derivatives at n = 1 (for all k ^ /cft)- 
Then Eq. ( ^Q]) implies that hp^{k, n, m) has the same property as a function of m at m = 0. 
Hence the expression ( ^0]) may be used on the right-hand side of Eq. ( |U1| ) for small positive 
m, up to an error of O(m^). Inserting the appropriate polynomials for momenta in the 
interval (n + 3m) /4 < k < min(3(n — m)/4, 1 — (3n -f m)/4) we obtain 

4?2(fc,H = E (^^!)(i+"^r'^^-^'a.+2(^ (C2) 

r=0 V / 

which holds for all k and small m. 

The following definitions will help keep the notation compact. Let P be any of A^^\ A^'^\ 
^(2)^ ^(3)_ jj-^ addition to the polynomials Pp{x, y), we define Pp{x, y) = Pp{x, y) + Cp_i(l, y), 
and introduce their generating functions 



P{x,y,z) = Y,Pp{^^y){-^y^ (C3) 

p=2 

P(") (x, z) = P{x, y, z) - illzll + 1 In 1^11 . (C4) 



2 2 \-z 



Now we are ready to collect the relations between Eqs. (|B^ ) and ( p5|) and between ( pSf) 
and (p^ ) that Eq. ( |6l]) provides [see Appendix Q. Setting x = k/n, y = m/n, and s = 1/n 
in these relation, we obtain (for arbitrary s) 

A%{x, y) = Y. (^) - yy^\-sr-'-A%{^^, fEi), (C5a) 

r=0 ^ ^ 

^'il^, 2/) + 4?2(^ - 2/) = E (^) - 2/)^^'(-^)"~^'^?+2(l^' 7^1 (C5b) 

P 



i?g2(^, + 4?2(^ - ^, = E (^) - l/)'''''(-^)"""4?2(i^, fE^), (C5c) 



BpUi^^y) + Bp}2{s -x,y) = -Cp+i{l,y) -sCp{l,y), (C5d) 



and other relations that are in fact implied by these. By repeated combination of Eqs. (|C5|) 
we find the important equality 

A« {x, y) + Af^ (x, y) = A^ {x,y) + 5f {x,y). (C6) 
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Furthermore, by setting s = 2x in Eq. ( |C5c| ) we immediately obtain the exphcit expression 

Bpl2{x,y) = -Cp+i{l,y) - xcp{l,y). (C7) 
In terms of generating functions our results so far can be expressed as 

A^'\x,y,z) = Qo{^,^)+O{y'), (C8) 

I(^) (x, y,z) = il~ sz) 1(1) fEj, ^) , (C9) 

1(1) (x, z) + A^'^ {s - X, ^) = (1 - sz) 1(2) ^) , (CIO) 

5(2)(x, z) = (1 - sz) 5(2) ffi, ^) = f In (Cll) 

1(3) (x, z) = 5(2) (x, + (1 - .z) 1(3) (fzf , frj, ^) , (C12) 

We proceed to determine A'^'^\ First we obtain y4(3)(x, 0, 2;) from Eq. (|C10|) at ?/ = by 
inserting the expansion ( p8|) and differentiating with respect to s; we then set s = 1, so that 
the error term vanishes. This yields an expression for (iy4(3)(x, 0, 2)/(ix, which we integrate 
with respect to x, using A(3)(3/4,0,z) = 0. This result for A(3)(x,0, z) is used in Eq. ( P12|) 
with s = 1, which yields A^^\x^ y, z) for arbitrary 



^(3) ( ) ^ {l-yz){Ax-l-?>y) . 



■Qi 



4 V l-i/ ' i-yz 

2(1 - 2/2;) 



The dot indicates derivative with respect to second argument. Next we derive A^^\ First 
we obtain A^^\x, 0, z) from Eq. ( |C10D at s = 1 and y = 0, using (|C^ ) at m = 0. This result 
for A(i)(x,0, 2;) is used in Eq. ( |C9| ) with s = 1, which thus yields A^^\x,y, z) for arbitrary 

y, 

A^'\x^y^z) = (i^l^H^Jl^f^Q^fl^, (iz^ 

, (l~3?/)2: + 2 / 43:-2+2y (l-y)z ' 
/^4£-2+2y 

2(1 -yz) 'V ^-y ' 1-?^- 

+ i?(2)(x,y,z). (C14) 

Finally, A^'^\x,y,z) is obtained from Eq. (|C6|). These results can be rearranged into 
Eqs. dD-dl). 
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FIG. 1: Double occupation d{g,n,m) as a function of the variational parameter g for various 
densities n and magnetizations m. Note that d{g,n,m) = d{g,n,0) — d{g,m,0). 
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FIG. 2: Momentum distribution Uka of (a) majority and (b) minority spin electrons, as a function 
of A; for 5 = 0.1 and density n = 1. 
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FIG. 3: Momentum distribution Uka of (a) majority and (b) minority spin electrons as a function 
of A; for 5 = 0.1 and density n = 0.8. 
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FIG. 4: Discontinuity of the momentum distribution at the Fermi vector for densities (a) n 
1 (in this case = q^) and (b) n = 0.8. 
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FIG. 7: Critical interaction Uc{n) for the transition from a paramagnetic to a fully polarized 
ferromagnetic state for the Hubbard chain with nearest-neighbor hopping as obtained from the 
Gutzwiller wave function. The Hartree-Fock result, [Eq. (|73D], is also shown; it has a shallow 
maximum at n ~ 0.856. Note that Uc{n) = Uc{2 — n). 
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FIG. 8: A class of Fcynman diagrams contributing to h^^{k,n,m). Solid (broken) lines repre- 
sent majority (minority) spins and vanish unless they carry momentum in the interval [— fepj, fcpt] 
([— fcpj., /cpj,]); see MV for diagrammatic rules. For k > S/cpj, the diagram in (a) vanishes since 
A; + A;i + A;2 G [— fepj, kpi] cannot be fulfilled. For k > Sfcpj, the diagram in (b) vanishes for similar 
reasons, and so on for higher odd multiples of fcpj^. 
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